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Abstract. We investigate the cosmological behavior of mimetic F(R) gravity. This scenario 
is the F(R) extension of usual mimetic gravity classes, which are based on re-parametrizations 
of the metric using new, but not propagating, degrees of freedom, that can lead to a wider 
family of solutions. Performing a detailed dynamical analysis for exponential, power-law, 
and arbitrary F(R) forms, we extracted the corresponding critical points. Interestingly 
enough, we found that although the new features of mimetic F(R) gravity can affect the 
universe evolution at early and intermediate times, at late times they will not have any 
effect, and the universe will result at stable states that coincide with those of usual F(R) 
gravity. However, this feature holds for the late-time background evolution only. On the 
contrary, the behavior of the perturbations is expected to be different since the new term 
contributes to the perturbations even if it does not contribute at the background level. 
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1 Introduction 

In order to explain the late-time universe acceleration one can follow two main directions. 
The first is to introduce the concept of dark energy in the framework of General Relativity 
(for reviews see [1, 2]), while the second is to modify the gravitational sector itself (for 
reviews see [3, 4]). The latter approach has an additional motivation, namely to improve the 
UltraViolet behavior that arises from the non-renormalizability of General Relativity and the 
difficulties of its quantization [5]. However, we mention that one can transform between the 
above directions, partially or completely, or construct various combined scenarios such are 
those with nonminimal couplings [6]. 

In order to construct gravitational modifications one usually adds higher-order correc¬ 
tions to the Einstein-Hilbert action. Amongst them the simplest model is that of F(R) 
gravity, where one replaces the Ricci scalar R in the action by an arbitrary function F(R) 
[7], which proves to have interesting cosmological implication such is the successful descrip¬ 
tion of inflation [8] (see [9, 10] for the analysis of the cosmological density perturbations), 
of late-time acceleration [11-32], or of both in a unified picture [33-35]. Furthermore, other 
higher-curvature models are those using the Gauss-Bonnet term G [36, 37] or functions 
of it [37-39], Lovelock combinations [40, 41], Weyl combinations [42, 43], Galileon mod¬ 
ifications [44-47], higher spatial-derivatives as in Horava-Lifshitz gravity [48-64], suitable 
self-interacting gravitational terms as in nonlinear massive gravity [65-68] etc. 
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One interesting class of gravitational modification, is that of mimetic gravity [69-76]. In 
these constructions one parametrizes the metric using new, but not propagating, degrees of 
freedom, and thus he obtains modified field equations which may admit a wider family of so¬ 
lutions. Usually, one can obtain solutions with an extra term proportional to a -3 , and that is 
why many authors talk about “mimetic dark matter”, i.e. a matter-like term of gravitational 
origin. In these lines, in [77] the authors added an F(R) modification in the framework of 
mimetic gravity, and showed that the resulting cosmology can accept new solutions comparing 
to usual F(R ) gravity or usual mimetic gravity. Thus, “mimetic F(R) gravity” corresponds 
to a new class of gravitational modification that deserves further investigation. 

In the present work we are interested in studying in detail the cosmological behavior 
is scenarios governed by mimetic F(R) gravity. In order to bypass the complexity of the 
involved equations that do not allow for analytical solutions, we use the powerful method 
of dynamical analysis, which extracts information about the global behavior of the scenario 
[78, 79]. However, a significant difference comparing to usual mimetic gravity is that now 
in the Friedmann equations, apart from the term proportional to a -3 , we obtain the F(R) 
contributions. Since both these contributions can only be observed through gravitational 
observations, it is impossible to separate them, and hence one must include them in a unified, 
dark-energy sector. The situation is similar to the case of “dark radiation”, i.e. a term of 
gravitational origin proportional to a -4 , that appears in many models, which is usually 
considered a part of the dark energy sector [80-84] , even if in this case one can use Big Bang 
Nucleosynthesis in order to constrain it. Therefore, in the scenario at hand, where such a 
constrain is moreover absent, the incorporation of the new terms in an effective dark energy 
sector is the only self-consistent approach. 

The plan of the work is the following: In Section 2 we review the scenario of mimetic 
F(R) gravity and we apply it in a cosmological framework. In Section 3 we investigate the 
dynamics in the case of an exponential F(R) form, while in Section 4 we perform the analysis 
for a power-law ansatz. In Section 5, for completeness, we provide the tools for a general 
analysis for arbitrary F(R) forms. Finally, in Section 6 we discuss the physical features of 
the obtained results, while section 7 is devoted to the Conclusions. 

2 Mimetic F(R) gravity and cosmology 

In this section we provide a brief review of mimetic F(R) gravity [77]. As we mentioned 
in the Introduction, the idea behind the general class of mimetic gravities [69-76] is that 
parametrizing the metric using new (but not propagating) degrees of freedom one can obtain 
modified field equations which may admit a wider family of solutions. For instance, after the 
action of a metric theory is given, a convenient parametrization of the metric g IJL/ is [69] 

9fiu = -g^dp^d^g^ , ( 2 . 1 ) 

and thus the action variation will be performed in terms of both g /JU and cj) (for an equivalent 
formulation using Lagrange multipliers see [70, 71, 85]). We stress that relation (2.1) implies 
that 

g(g f iuAT u d ll 4>du(f > = - i, ( 2 - 2 ) 

which shows that the scalar field will not be a propagating degree of freedom [69, 77, 85, 86]. 
Additionally, due to the above parametrization, the mimetic extension of the initial theory 
has become conformally invariant. In summary, variation with respect to g will give rise to 
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the traceless part of the Einstein equations, while variation with respect to <j> gives the trace 
part of Einstein equations modified by an extra prefactor and thus allowing for a wider class 
of solutions. Such solutions may have an effective dark-nratter-like component, and since in 
some sense the whole theory mimics a dark matter sector, the theory is named “mimetic” 
gravity. Lastly, note that g^ does not appear in the final equations of motion, since it can 
be eliminated in terms of the initial metric g, JU and (j). 

Let us apply the above general instructions in the usual metric F(R) gravity following 
[77]. We start from the standard F(i?)-action 


S = 



' F(R) 

2 K 2 


+ F r 


(2.3) 


where k 2 is the gravitational constant, R is the Ricci scalar calculated by the metric and 
C m stands for the matter Lagrangian. Parametrizing the metric according to (2.1) we obtain 


S = j d A x\J—g (g^,(j>) 


F{R{g,xvA)) 


2k 2 


+ Fr, 


Hence, variation with respect to g lw gives [77] : 


(2.4) 


2 9^yF ( R (< ipuvi ( t > )) R (g^, (f>)F r ( R (g^vi </*)) 


+v {g,iv; <t>)^ v (g (g^, Fr (R (g, w , </>)) 


-g ( gnu, 4>) uv □ GW, 4>) f r {R (<W> 4>)) + 


Y)n U 11 Y) ± H YJJ T ™ 

+d ll (j)d v (j) 2 F (R (<W> </>)) - R (<W> 0) Fr (R (<W, 4>)) 

-2r 


-3D [g(g^u, 4*)^ ) Fr (R (<W, 4>)) + k 2 T 


= 0 , 


(2.5) 


where Fr stands for dF(R)/dR, and □ are respectively the covariant derivative and box 
operators with respect to g (first expressed in terms of g flll and its derivatives and then 
expanded in terms of g^,^ through (2.1)), and T^ u is the matter energy-momentum tensor 
arising from C m . Additionally, variation of (2.4) with respect to (j) leads to 


V (g (WWV^ \py,4> 2 F(R(g^,4))) - R (g jlv , (f>) Fr (R (<W, <f>)) 
-3D (g {guv, <j>)^ Fr (R (<W, 4>)) + k 2 T } = 0 . 


( 2 . 6 ) 


with T = g {gnu- (p) lw T /IU the trace of the matter energy-momentum tensor T^ u . Since the 
above equations do not contain g^ u explicitly, but only through the combination (2.1), in the 
following we omit the g and <f> dependence of the various quantities. 

We mention that every solution of standard F(R) gravity is a solution of the above 
mimetic F(R) gravity, however the opposite is obviously not true. Apart from the wider 
class of solutions, the advantage of the new theory is that it is conformally invariant. 

Since we are interested in investigating the cosmological implications of mimetic F(R) 
gravity, in the following we consider the flat Friedmann- Robertson-Walker (FRW) metric 


ds 2 = — dt 2 + a 2 (t) 8ijdx l dx\ 


(2.7) 
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where a(t) is the scale factor. Since (j> is homogeneous in this case, the constraint (2.2) 
leads to <j) = t, which simplifies significantly the equations. Moreover, the Ricci scalar as 
usual becomes R = 6 H + 12 H 2 , where H = a/a is the Hubble parameter and dots denoting 
differentiation with respect to the cosmic time t. Under the above considerations, the 00 and 
ii components of (2.5) lead to the same equation 

0 = 2 F rrr R 2 + 2 F rr R + 4 HF rr R -2 (h + 3 H 2 ^j F R + F(R) + 2 K 2 p m , (2.8) 


while (2.6) gives 

= 2 F(R) — 6 (h + 2H 2 j F R + 3 F rrr R 2 + 3F rr R + 9 HF rr R + k 2 (3 p m — p m ) , (2.9) 

where C is a constant of integration, and p m and p m are respectively the energy density 
and pressure of the perfect-fluid matter sector (T = — p m + 3 p m ). These are the Friedmann 
equations of the scenario at hand. They can be rewritten as 


C<j> _ k 2 Pm F{R) HRFrr 

3 a 3 F R 3 F r 6 F R F R 

C HF RR RFrr 2 H 2 F r F{R) n 2 (3p m + p m ) 

3 a 3 R 2 R R 2 R 2 3 R 2 3 R 2 


(2.10a) 

(2.10b) 


where we have reduced the Raychaudhuri equation (2.8) to its simpler form (2.10a) by elimi¬ 
nating the third-order derivative F RRR through (2.9). Finally, for the purpose of the following 
analysis, it proves convenient to re-express the first Friedmann equation as 


H + 


RFrr 
2 F r 


F(R ) 
6 F r 


R 2 Pm C# R R 2 F rr 2 
3 F r + 3 a 3 F R + 6 + 4 F R 2 


( 2 . 11 ) 


In the case where = 0 we re-obtain the equations of motion of standard F(R) gravity, 
however in the general case we obtain a correction-term proportional to a~ 3 . Hence, once 
again we verify that all solutions of standard F(R) gravity are solutions of the above theory 
for Cff, = 0. 

The extra term proportional to a -3 is present in all mimetic gravity versions [69-77] 
and since it mimics an effective matter sector it gave the name “mimetic” to this class of 
theories. However, since this term is an effective term of gravitational origin, and thus it 
will not appear in the future experimental (direct) verification of dark matter, and since in 
the Friedmann equations in appears alongside the F(R) terms, in the present work, and in 
contrast with the usual mimetic considerations, we prefer to incorporate it inside the effective 
dark energy sector. The situation is similar to the case of “dark radiation”, i.e. a term of 
gravitational origin proportional to a~ 4 that appears in many models. Dark radiation is 
considered a part of the effective dark energy sector and not a part of the radiation sector 
(the physics of which is more or less known), although one can use Big Bang Nucleosynthesis 
data in order to constrain it independently of the rest dark energy sector [80]. Hence, in the 
present model we will consider the term U^/a 3 as part of the effective dark energy sector, 
alongside with the F(R) terms, although one might use gravitational lenses data in order to 
constraint it independently of the F(R) terms. 
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In these lines, we can rewrite the Friedmann equations (2.8),(2.9) in the usual form 


H 2 = —(p m + PDE ) (2.12a) 

k 2 

H = ——{pm+Vm + Pde + Pde) , (2.12b) 


defining the energy density and pressure of the effective dark energy sector as 


Pde 


1 

K 2 


RFr-F 

2 


3 HRFrr + 3 H 2 (1 - Fr) + % 


(2.13a) 


Pde 


R 2 Frrr + 2 HRFrr + RFrr + 


F - RFr 


- (2H + 3H 2 ) (1 -Fr) 


. (2.13b) 


Hence, as we can see, the effect of the mimetic F(R) gravity comparing to the standard 
F(R) one, is an extra term in the effective dark energy density, while its pressure remains 
unaffected. Additionally, we can define the dark energy equation-of-state parameter as usual 
as 


Wre 


Pde 
Pde ’ 


(2.14) 


while the corresponding quantity of the matter sector is w m = Pm/Pm, which satisfies — 1 < 
w m <• 1- We mention here that we wrote the above Friedmann equations and we defined 
the dark energy quantities using the initial gravitational constant k 2 and not the effective 
one k 2 /Fr, in order to ensure the separate conservation of the dark energy and the matter 
sectors [7], namely 


Pde + 3F[(pde + Pde) = 0, (2.15a) 

Pm + 3H(p m +Pm) = 0. (2.15b) 

Lastly, as usual, the F(R) form is forced to satisfy the following general conditions [87]: 
The existence of a stable Newtonian limit requires 


| F(R) — i?| <C R, I Fr — 1| •C 1, RFrr <C 1, (2-16) 


in order for the non-GR corrections to a space-time metric to remain small (the last condition 
implies that the Compton wavelength is much less than the radius of curvature of the back¬ 
ground space-time) [87]. The ghost avoidance and classical and quantum stability requires 
[88, 89] (see also [90]) 

Fr > 0, Frr > 0. (2-17) 


Note that if Frr becomes zero for a finite R = R c , then a weak (sudden) curvature singularity 
is generally formed [87]. In the absence of matter, the asymptotic future stability of the de 
Sitter solutions requires 


Fr\r=r 1 

Frr\r=r 1 


> Ri, 


where R\ satisfies RFr — 2 F(R) = 0 [91]. 

Similarly to the standard F(R) case, one can in principle impose the desired a(t) be¬ 
havior and suitably reconstruct the F(R) form that generates it [77]. However, in this work 
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we are interested in the inverse procedure, that is first consider a specific F(R) form and 
then investigate the induced universe evolution. In order to achieve this independently of the 
specific initial conditions, in the following section we apply the powerful method of dynam¬ 
ical analysis [92—99] . In particular, we first transform the cosmological equations into their 
autonomous form and we extract the corresponding critical points. Then, we linearize the 
perturbations around these critical points, and we express them in terms of the perturbation 
matrix. Hence, the eigenvalues of this perturbation matrix for each critical point, determine 
its type and stability. 


3 Mimetic F(R) gravity with exponential form 


In this section we examine the behavior of mimetic F(R) gravity, under an exponential F(R) 
ansatz of the form 

F(R) = A [exp ( pR ) — 1] , (3-1) 

which is one of the most well-studied in standard F(R) gravity [7]. For convenience we focus 
on the physically interesting A > 0 and p > 0 cases, although the investigation of the general 
case is straightforward. Additionally, we parametrize this F(R) form as F(R) = f(R) — A, 
with f(R) = A exp ( pR ), and we choose A = p~ l . 1 

The Friedmann equations can now be expressed as: 


2 = ^Prn _J__ HRfRR R A 

3 a 3 f R 3 f R 6 f R f R 6 6f R ' 

H = -H 2 - _A_ - ^ + _L + _ A 

3 a 3 f R 3 f R 6 f R f R 6 f R 

while equation (2.10b) becomes 


(3.2a) 

(3.2b) 


R = -Stsg + 2 FL _ Hk _ ^ (^ R -1) _ ^ (3 . 3) 


3 pa A 


P 


3 p 2 


3 p 


In order to transform these equations into their autonomous form, we need to introduce 
suitably defined auxiliary variables [92-99]. Thus, we define the normalized variables 


p _ C 4) H _Rf RR f _ n 2 Pm 

3(rFD2 /fl ’ ^ D ' x 2Df R ’ !J 6D 2 f R ’ Z ZD 2 f R ' 


where 


D = 


\ 


H + 


+ 


/ 


R/rr \ 

2/k J ' 6 f R 


= \I(H + - P R ] +-. 


(3.4) 


(3.5) 


Moreover, we define two more auxiliary variables, which in the present example are related, 
namely 


- Rfri p 

r — -— = —pR, 


m 


J 

RJrr 

fn 


= pR = —r. 


(3.6) 


1 Note that at late times, i.e. for small curvatures (R <fC 1), we have F(R) ~ R + ^R 2 + A + 0(R 3 ). 
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The role of these two variables will become clear in section 5. 

From the definitions (3.4) we immediately extract the constraint (Q + x ) 2 + y = 1, while 
the first Friedmann equation (3.2a) leads to the constraint P + x 2 + y(e r — r) + z = 1. In 
order to handle this transcendental expression we introduce the additional variable 


Ha 


A 


QD 2 f R 


= e y = e 


i — (Q + %Y 


i.e., we have the restriction 


= 0 . 


Ha - e r 1 — (Q + x ) 2 
Additionally, the constraint equation becomes 

P + x 2 — ry + + z = 1. 


(3.7) 


(3.8) 


(3.9) 


The appearance of the above two constraints and the constraint (Q + x) 2 + y = 1 allows 
us to eliminate three auxiliary variables, for instance y and z and r. Hence, the dynamical 
equations for the remaining variables write as 


P' = P {x - (Q + x) {3r [(Q + x) 2 — l] + 2 Qx — x 2 + 3Ha} } 

- 3 Pw m (Q + x) [P + r[(Q + x) 2 - l] + x 2 + Ha - 1} , (3.10a) 


Q — — Q3(9 r + 2)x + Q~ [r (9ar — 5) + x 2 + 3Ha + l] 

— Qx [r (7 — 3x 2 ) + x 2 — 3Ha + 3] —2 r ( x 2 — l) j 
— ‘^QwmiQ + x )|P + r [(Q + x) 2 - l] + x 2 + Ha — 1 j, (3.10b) 

x' = — ^{3 Q^ rx — Q 2 [—(9r + 2)x 2 + r + 4] + Qx [r (9x 2 — 5) + x 2 + 3Ha — 5] 

+ ( x 2 — l) [r (3x 2 — l) — x 2 + 3Ha — 3] | 

- ^w m [x(Q + x) - 1] |P + r [(Q + x) 2 - l] + x 2 + Ha - l|, (3.10c) 

Ha = — Ha j3Q 3 r + Q 2 (9r + 2)x + Q [r (9x 2 — 3) + x 2 + 3 Ha — 3] 

+ x [3r (x 2 — l) — x 2 + 3 Ha — l] | 

— 3Ha Wm(Q + 2: )^-f > + r \(Q T x) 2 — l] + x 2 + Ha — 1 j - , 
where r is expressed in terms of the other variables as 


(3.10d) 

(3.11) 
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In the above equations, the primes denote derivatives with respect to the new time variable 
ry defined as drj = Ddt. Thus, the system (3.10) determines a flow on the region of the phase 
space defined as 2 


'hi = 


Q,x, Ha) : |Q + x| < 1,0<.P + x 2 — r [ 1 — (Q + x) 2 ] + 12 a A 1; 


r = In 


r^A 


1 — (Q + x) z 


,p> o}. 


(3.12) 


Lastly, in terms of the auxiliary variables (3.4) and r, explicitly given by (3.11), the mat¬ 
ter and dark-energy density parameters from (2.12a),(2.13a), the deceleration parameter, the 
dark-energy equation-of-state parameter (2.14), and the total equation-of-state parameter, 
are written as 


K 2 p m _e r {l-P + (e r -r) [(Q + x) 2 - l] - x 2 } 

' 'm. — 


m 3H 2 

_ k 2 p DE e _r Ai 

= 


Or 


3 H 2 Q 2 ’ 

1 H _ r[l-(Q + x) 2 ] 

q = - l ~w - (p +1 - 

e r {Q 2 (2r + 3w m - 1) + Q(4rx + 6 xw m ) + (x 2 - l) (2 r + 3 w m )} 


wde = 


w r 


3Ai 

{P + r [(Q + x) 2 - l] + x 2 - 1} 


Wtot = -J- - 


i 2 H _ 2q - 1 _ 2r [(Q + x) 2 - l] 1 

-1 o tto ^ r 77, 


3 H 2 3 


3 Q 2 


(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 


where Ai = e r (2Qx + x 2 — l) — {P + r [(Q + x) 2 — l] + x 2 — l}. 

The scenario of mimetic F(R) gravity with the exponential form (3.1), i.e. the system 
(3.10) that lies on the reduced phase space (3.12), admits three isolated physical critical 
points (note that the appearance of the constraint (3.11) reduces their number significantly), 
which are displayed in Table 1 along with their existence and stability conditions. The details 
of the analysis and the calculation of the various eigenvalues of the 5x5 perturbation matrix 
are presented in Appendix A. Furthermore, for each critical point we calculate the values of 
various observables, such as the density parameters, the deceleration parameter and the dark- 
energy and total equation-of-state parameters, given by (3.13)-(3.17), and we summarize the 
results in Table 2. 


4 Mimetic F(R) gravity with power-law form 

In this section we study the behavior of mimetic F(R) gravity under a power-law F(R) ansatz 
of the form 

F(R) = R + aR n - A, (4.1) 

which is also one of the most well-studied in standard F(R) gravity [7]. We focus on the 
physically interesting A > 0 case, although the analysis of the general case is straightforward. 
Furthermore, we parametrize this F(R) form as F(R) = f(R) — A with f(R) = R + aR n . 

2 Note that from (3.10a) it follows that the sign of P is invariant, and recall that we have assumed that 
fn> 0 which implies that P > 0. 
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Name 

P 

Q 

r 

X 

Oa 

Existence 

Stability 

Si 

0 

0 

0 

0 

1 

always 

nonhyperbolic (see numerics) 

s 2 

0 

Qci 

2Q^ 

A 

0 

l - 2 Q 2 cl 

always 

saddle 

S 3 

0 

Qc2 

2Q^ 2 

*3c2 —3 

0 

1 - 2 Q 2 c2 

always 

saddle 


Table 1. The real critical points of the system (3.10) of mimetic F(R ) gravity with the exponential 

form (3.1) and their existence and stability conditions. The parameters Q c \ and Q c2 correspond to 

2 qg 

the two roots of the transcendental equation 2 Q\ — e Q = -1 (Q 2 — l) — 1 = 0, which numerically are 
found to be Q c \ = —Q c 2 « 0.666, which belong to the interval ( — ^ j and thus > 0. 


Name 


^DE 

q 

WDE 

Wtot 

Si 

arbitrary 

arbitrary 

arbitrary 

arbitrary 

arbitrary 

s 2 

0 

1 

-1 

-1 

-1 

S3 

0 

1 

-1 

-1 

-1 


Table 2. The real critical points of the system (3.10) of mimetic F(R) gravity with the exponential 
form (3.1), and the corresponding values of the matter and dark energy density parameters, of the 
deceleration parameter, and of the dark-energy and total equation-of-state parameters, calculated 
through (3.13)-(3.17). 


The Friedmann equations can now be expressed as: 


2 = _ J_ _ HRf RR R A 

3 a 3 f R 3 f R 6 f R f R 6 6 f R 

H = -H 2 - -A_ - ^ + A + _ A 

3 a 3 f R 3 f R 6 f R f R 6 f R 

while equation (2.10b) becomes 



c<t>R 2 - n , U 2 ' 

3a(n — l)na 3 
k 2 (3 w m + 1 )p m R 2 ~ n 
3 a(n — l)n 


2 R 2 ~ n 2 R 

a{n — 1 )n n — 1 


- HR - 


A R 2 ~ n R 3 ~ n 

3 a(n — l)n 3a(n — 1 )n 


(n - 2 )R 2 
R 
R 2 

3 (n — l)n" 


(4.2a) 

(4.2b) 


(4.3) 


In order to transform these equations into their autonomous form we introduce the normalized 
variables 


with 


P = 


C<p 

3a 3 D 2 f R 


Q 


H _ Rf RR _ f _ K 2 p m 

D 1 2 “ 2 Df R ' V ~ 6D 2 f R ’ ' " 3D 2 f R ’ 


D 


N 


H + 


R/rr 
2 fR 


2 

+ A 
6 //? 


(4.4) 


(4.5) 
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Moreover, we define two additional auxiliary variables, which in the present example are 
related, namely 


Rf R R ( anR n 1 + l) 

/ “ aR n + R 

R/rr n(l + r) 

m = — — = —-. 

JR r 

(4.6) 

Finally, similarly to the previous section, we define 


6 D^f R 

(4.7) 

Hence, from the definitions (4.4) and the first Friedmann equation (4.2a) 
above auxiliary variables satisfy the constraints 

we deduce that the 

P + x 2 -ry + z + fl a = 1, 

(4.8) 

and 


(Q + x) 2 + y= 1. 

(4.9) 


Using the above two constraint equations in order to eliminate two auxiliary variables, namely 
y and z, we finally result to the following autonomous dynamical system: 


Px 3 (3 nr 2 + 2r 2 + 4nr + n) + PQx 2 [4r 2 + n(r + l)(9r + 5)] 

n(r + 1) 

Px [2r 2 + 3n(r + l)(r — Ua + 1)] — PQ 2 x [2r 2 + n(r + l)(9r + 4)] 

n(r + 1) 

+ 3 PQ [(r - Ua) - Q 2 r ] - w m ^3Px 3 (r + 1) - Px [3(r - f^y + 1) - 9 Q 2 r - 3 P] 

+ PQx 2 (9r + 3) — 3 PQ (r — Ha + 1 — Q 2 r — Q) |, (4.10a) 

^ Qx 3 [(3n + 2)r 2 + 4nr + n] + Q 2 x 2 [4r 2 + n(r + l)(9r + 5)] 

^ 2n(r + 1) 

Qx [2r 2 + n(r + l)(7r — 3 Ua + 5)] — Q 3 x [2r 2 + n(r + l)(9r + 4)] 

+ 2 n(r + 1) 

3rQ 4 Q 2 n 2 

---1- — (5r - 3S2 a - 1) + rx -r 

1 ) 

-{(r + 1 )x 3 Q + (3r + l)x 2 Q 2 + x [3 rQ 3 + PQ - Q(r - Ua + 1)] 

+rQ 4 + PQ 2 -Q 2 (r-n A + 1)} , (4.10b) 


, 2 r(n + r)x 

n 


(4.10c) 


10 














x 4 [(3n + 2)r * 2 + 4nr + n] + Qx 3 [4r 2 + n(r + l)(9r + 5 )] 

2 n(r + 1) 

x 2 { [2r 2 + n(r + l)(4r — 312 a + 4)] — Q 2 [2r 2 + n(r + l)(9r + 4)] } 

+ 2n(r + 1) 

+ [(5x — 30 a + 5) — 3Q^r] + ^-(r + 4) — — (r — 30a + 3) 

Qi ij 

4- {(r + l)x 4 - Q(3r + l)x 3 + x 2 [-3rQ 2 - P + 2(r + 1) - 0 a] 

-Qx ( rQ 2 4- -P - 3r 4- Oa - l) 4- P + Q 2 r -r + Oa — l} , (4.10d) 

Oax 2 { [4r 2 4- n(r 4- l)(9r 4- 5)] Q + [(3n 4- 2)r 2 4- 4nr 4- n] x} 


Oax { [2r 2 4- 3 n(r + l)(r — Oa 4- 1)] — Q 2 [2r 2 4- n(r 4- l)(9r 4- 4)] } 

n(r + 1) 

— 30aQ ( rQ 2 — r + Oa — l) — 30 a w m {x (3rQ 2 4- POa — r 4- Oa — l) 

4-(r 4- l)x 3 4- Q \rQ 2 4- (3r 4- l)x 2 + P — t4-0a— l] } . (4.10e) 

In the above equations the primes denote derivatives with respect the new time variable r/ 

defined as d?? = Ddt. Hence, the system (4.10) defines a flow on the region of the phase space 

3 


^2 := {(P 1 , Q,r,x, Oa) :0<P-r[l-(Q4 x) 2 ] 4- x 2 4- Oa < 1, \Q + x| < 1, P > 0} . 

(4.11) 

Finally, the matter and dark-energy density parameters from (2.12a),(2.13a), the de¬ 
celeration parameter, the dark-energy equation-of-state parameter (2.14), and the total 
equation-of-state parameter read as 


_ K 2 p m _{n- 1 )r {P + r [(Q + x) 2 - l] + x 2 + 0 A - l} 

' “m. — 


0 de — 


3 H 2 
2 

K PDE 


Q 2 (n + r ) 


Ao 


3 H 2 Q 2 (n + r ) 

= = r[l-(Q + x) 2 ] i 


WDE = 


H 2 Q 2 + 

(n 4- r) [Q 2 (2r — 1) 4- 4 Qrx + 2 r (x 2 — l)] 


3A, 


+- 


w 7 


i(n - 1 )r {P + r [(Q 4- x) 2 - l] + x 2 + Oa - 1} 


A 2 


W t ot = “I - 


2 H _ 2q - 1 _ 2r [(Q + x) 2 - l] ( 1 

1 o TTO o o7v? f" TT 1 


3 H 2 


3 Q 2 


(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 


where A 2 = (n — l)r [P 4- 2 Qrx + (r 4- 1) (x 2 — l) 4- Oa] 4- Q 2 (r 4- 1) [ n(r — 1) — r\. 

The scenario of mimetic F(R) gravity with the power-law form (4.1), i.e. the system 
(4.10), admits 14 x 2 4- 1 = 29 isolated physical critical points and three curves of critical 

3 Note that from (4.10a) it follows that the sign of P in invariant, and recall that we have assumed that 
fa > 0 which implies that P > 0. 
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Name 

P 

Q 

r 

X 

n A 

Existence 

Stability 

T{ 

0 

0 

0 

€ 

0 

always 

unstable (stable) 

n 

0 

0 

—n 

€ 

0 

always 

saddle 

T 3 

0 

2e 

0 

—e 

0 

always 

saddle 

n 

0 

2e 

—n 

—e 

0 

always 

vmstable (stable) for w m < | and n < 1 or n > | 
saddle otherwise 

n 

9 

2 e 

3 

0 

3 

0 

always 

saddle 

Ti 

9 

ze 

3 

—n 

3 

0 

always 

saddle 


Ti 

2ns 

y/ n(n+2)+3 

—n 

3(n-l) £ 

y/ n(n+2)+3 

0 

| < n < ± (i3 + yra) SJ 1.35 

Non-stable 

ti 

0 

(In - l)r 2 

—n 

(n - 2)r 2 

0 

|<n<lorn>| 

stable (unstable) for w m > — l,n > 2 
saddle otherwise 


0 

2e 

0 

(W-l)e 

0 

-1 < < § 

saddle 

Tb 

0 

2e 

—n 


0 


saddle 

n 1 

0 

2 n 

VTZ 

—n 

%/F 3 

0 

5 < n < 2 1 < Wm < ~ 8n f+ 13n ~ 3 

4 n Z, 1 ^ w m A 6 n a - 9 n +3 

or 0 < n < -1 < w m < §(4n — 3) 

or n — 2, w m = — 1 

stable (unstable) for 
l<n<2,-l<w m < 

or n — 2, Wm — — 1 
saddle otherwise 

T{ 2 

0 

vie 

2 e 

-2 

0 

0 

always 

non-hyperbolic with 4D stable (unstable) manifold for 

0 < n < 2, — 1 < w m 

n 3 

0 

\J 

—n 

0 

2—n 
n+2 

n > — 2 

numerical determination (see Appendix B) 

T u 

0 

1 

r cl4 

0 

0 


saddle 

TlS 

0 

0 

0 

0 

1 

—1 < w m < 1 

non-hyperbolic 

Tie 

0 

Qc16 

2 

0 

Qc16 

1-C&6 

-3 < 0 cl6 < § 

saddle 

T 17 

0 

Qc17 

mn 

0 

1 — 2Q%17 

Ql 17 < 1 

stable for n < 0, — 1 < w m < 1, 5^ 75t ^ 32 < Q c n < 

or 0 < n < |, — 1 < w m < 1, Q c i 7 > -^= 
or n > §, -1 < w m < 1, ^ < Q c 17 < ^ 3^=2 

Tfs 

0 

3 

—n 

'2e 

3 

t> 

9 

always 

stable (unstable) for w m > — 1, 0 < n < 1 
saddle otherwise 


Table 3. The real critical points and curves of critical points of the system (4.10) of mimetic F(R) 
gravity with the power-law form (4.1). We use the notation e = ±1, where e — +1 corresponds to 
expanding universe and e = — 1 to contracting one, with the stability conditions outside parentheses 
corresponding to e = +1 while those inside parentheses to e = — 1. We have defined Ti = 2 6 ; 

y/n—le 


r 2 = 


and 1?3 = n 2 + 9(n — 1 ) 2 w^ + 6 \(n — 4)n + 2] w m + 2n + 3. 


yj n[n(9n—19) + 13] —4 

r c i 4 , Q c i 6 and Q c 17 are the parameters of the corresponding curves. 


n(n+2)+3 

Additionally, 


points (one of them, namely T 14 , exist only for a specific value of the parameter w m ), which 
are displayed in Table 3 along with their existence and stability conditions. The details of 
the analysis and the calculation of the various eigenvalues of the 5x5 perturbation matrix 
are presented in Appendix B. Furthermore, for each critical point we calculate the values of 
various observables, such as the density parameters, the deceleration parameter and the dark- 
energy and total equation-of-state parameters, given by (4.12)-(4.16), and we summarize the 
results in Table 4. Observe that for some specific points having either Q = 0 or r = —n the 
expressions (4.12)-(4.16) are not well defined (NWD), since the involved limits depend on 
the limit order. 

5 Dynamical analysis for general F(R) forms 

As we saw in the previous sections, in order to perform the stability analysis one needs to 
choose a specific F(R) ansatz. However, this is restricting since for different F(R) forms one 
must repeat the whole analysis from the start. Hence, in the present section, for completeness, 
we extend the usual procedure in order to be able to perform the analysis for arbitrary F(R) 
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Name 


^DE 

Q 

WDE 

w tot 

T{ 

NWD 

NWD 

NWD 

NWD 

NWD 

n 

NWD 

NWD 

NWD 

Wm 

NWD 

T e 

1 3 

0 

1 

1 

1 

3 

l 

3 

rpe 

1 4 

NWD 

NWD 

1 

NWD 

1 

3 

rpe 

1 5 

0 

1 

1 

i 

3 

1 

3 

npe 

1 6 

NWD 

NWD 

1 

NWD 

1 

3 

T e 

1 7 

NWD 

NWD 

—— 1 

2 n 1 

NWD 

~ l + n 

OOf^ 

NWD 

NWD 

1 + 3 1 
n—1 “ 1—2n 

NWD 

(7—6n)n+l 
6n' 2 — 9n+3 

rpe 

1 9 

0 

1 

1 

l 

3 

1 

3 

rpe 

1 10 

NWD 

NWD 

1 

Wm 

1 

3 

rpe 

1 11 

NWD 

NWD 

— 2n+3(Wm+l) 

2 n 

W m 

1 | ^m + 1 

' n 

rpe 

1 12 

NWD 

NWD 

NWD 

W m 

NWD 

Tf 3 

NWD 

NWD 

-1 

NWD 

-1 

Tu 

r c i4(l— n) 
n+r c ,14 

ra(r c i4+l) 
n+r c 14 

1 

i 

3 

1 

3 

T,5 

NWD 

NWD 

NWD 

NWD 

NWD 

Tie 

0 

1 

1 

l 

3 

l 

3 

Tn 

0 

1 

-1 

-1 

-1 

Tis 

NWD 

NWD 

1 

NWD 

1 

3 


Table 4. The real critical points and curves of critical points of the system (4.10) of mimetic F(R) 
gravity with the power-law form (4.1), and the corresponding values of the matter and dark energy 
density parameters, of the deceleration parameter, and of the dark-energy and total equation-of-state 
parameters, calculated through (4.12)-(4.16). We use the notation e = ±1, where e = +1 corresponds 
to expanding universe and e = — 1 to contracting one, with the stability conditions outside parentheses 
corresponding to e = +1 while those inside parentheses to e = — 1. NWD stands for “Not well-defined”. 


forms. Following the generalized method of [100], the idea is to suitably parametrize an 
arbitrary F(R) function and perform the dynamical analysis in general. Therefore, after 
this general analysis one can just substitute the specific F(R) form in the obtained results, 
without the need to repeat the whole dynamical elaboration from the beginning. 

In order to parametrize the arbitrary F(R) functions, we introduce the auxiliary vari¬ 
ables [16, 100] 

RFr 


RFrr 

Fr 


m = 


(5.1) 



































Furthermore, we introduce the normalization factor 


D = 


\ 


-id 


+ 


F 

6iV 


(5.2) 


and the normalized variables 


P = 


F 


n _ H_ _ RFrr _ _ 

3a 3 D 2 F R ’ I)' X ' 2 DFr V ~ 6 D 2 F R 


z = ./Lf™ , u = F r . (5.3) 


3 D 2 Fi 


R 


Since the consistency conditions require F(R) > 0 and F R > 0, it follows that P > 0 ,y > 
0 ,z > 0. Using the above auxiliary variables, the Friedmann equation (2.11) leads to the 
constraint 

P + x 2 — ry + z = 1, (5.4) 

while the definition of D gives rise to the additional constraint 


(Q + xf + y = 1. 


(5.5) 


Therefore, we can use the above two constraints in order to eliminate two variables, which 
for convenience are chosen to be y and 2 , through 


y = i- (Q + x ) 2 , 

z = 1 — P — x 2 + r\l — (Q + x) 2 ] 


(5.6) 


Defining a new time variable rj through dr] = Ddt, we can finally extract the autonomous 
form of the cosmological equations as 


p , = M(r) p 
r + 1 


— 4 Qx 2 + 2(1 — Q 2 )x — 2x 3 


Px (9 Q 2 r 2 + 11 Q 2 r + 4Q 2 — 3r 2 — 4r — 3) 


r + 1 


- 3rc m p|x [P + (3 Q 2 r - r - l)] + PQ + Q ( Q 2 r - r - l) + Q(3r + l)x 2 + (r + l)x 3 | 


PQ (9r 2 + lOr + 5) x 2 


r + 1 


+ 3P(1 — Q 2 )Qr — 


P (3 r 2 + 2r + 1) 
r + 1 


(5.7a) 


M (?’ 


Q' = - - -Q x (1 - x 2 ) - Qx (Q - 2x) 

r + 1 L 


+ 

3 

- w, 

2 


3 Q 4 r Q 3 [r(9r + 11) + 4]x 


2 2(r + l) 

{ Q 2 [P + r (3x 2 - l) + x 2 - l] + Qx [P + (r + 1) (x 2 - l)] + Q 4 + 3Q i rx) 


Q 2 f [r(9r + 10) + 5]x 2 
r + 1 

+ r (x 2 - l) , 


- 5r + 1 > + 


Qx {7 r 2 — [r(3r + 2) + l]x 2 + lOr + 5} 

2(r + 1) 


(5.7b) 
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X = 


— 2 Qx 3 + (1 — Q 2 )x 2 — x 
.2 


(9 Q 2 r 2 + 11 Q 2 r + 4 Q 2 — 4 r 2 — 6 r — 4) 


M(r) 
r + 1 
3 

- 2 w * 

+ (r + l)x 4 - P - Q 2 r + r + 1 - ^Qx (3 Q 2 r - 5r - 5) + ^ (Q 2 r + 4Q 2 - r - 3) 


x 


2 (r + l) 

(P + 3Q 2 r - 2 r - 2 ) + xQ (P + Q 2 r - 3r - l) + Q(3r + 1 ): 


Q (9 r 2 + lOr + 5) x 3 (3r 2 + 2r + l) x 4 

2(r + 1 ) 2 (r + 1 ) 

(5.7c) 

r' = 2 M(r)x, 

(5.7d) 

u = 2 xu. 

(5.7e) 

and we have the additional equation 



D' = 


M(r)Dx 


r + 
3 


— [ 2 Qx + (Q - 1)(Q + 1) + x 2 ] 


+ -Dw m [x (P + 3 Q 2 r - r - l) + PQ + Q ( Q 2 r - r - l) + Q(3r + l)x 2 + (r + l)x 3 ] 


+D 


x (9Q 2 r 2 + llQ 2 r + 4Q 2 — 3r 2 — 6r — 5) 3 / 9 N 

v - -— 7 -+ -Q (Q 2 r - r - 1) 

2(r +1) l ^ ’ 

Q (9r 2 + lOr + 5) x 2 | (3r 2 + 2r + l) x a 
+ 2(7+1) + 2(7+1) 


(5.8) 


where primes denoting derivatives with respect to rj, and with 

(5. 

assuming that m can be expressed as a function of r, namely m = m(r). Since the equation 
(5.8) is decoupled form the rest, we are allowed to investigate the restricted dynamical system 
defined in the phase space 

+ = {(P, Q, x, r, u) : \Q + x| < 1,0 < P — r [l — (Q + x) 2 ] + x 2 < 1, P > 0} . (5.10) 

Additionally, note that since the evolution equation for u is decoupled too, it follows that the 
Jacobian matrix of the extended dynamical system for (P, Q,r,x,u ) has an extra eigenvalue 
A u = jj^\x=x c ,u=u c i where (x c ,u c ) are the values of (x,u) at the equilibrium point. Hence, 
from (5.7e) we deduce the two limiting situations at an equilibrium point, namely [100]: 

• For x c = 0 it follows that \ u = 0. Thus, at the equilibrium point fn acquires a constant 
value, and the stability issue cannot be resorted by linear analysis. 

• For x c 0 it is required that u c = 0, which implies that fn = 0 at the equilibrium 
point. Additionally, X u = 2x c and thus perturbations along the rt-axis are conditionally 
stable in the extended phase space for x c < 0 . 

Lastly, the matter and dark-energy density parameters from (2.12a),(2.13a), the deceler¬ 
ation parameter, the dark-energy equation-of-state parameter (2.14), and the total equation- 


M(r) = 


r{ 1 + r + m) 


m 
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of-state parameter, can be expressed as 


_ K 2 Pm _ {P + r [(Q + x ) 2 — l] + x 2 - 1 } u 
= 3 h 2 ~ Q 2 ’ 


_ k 2 pde _ {p + r [( Q + x) 2 - l] + x 2 - 1} u + Q 2 
n DE = — 2 - - -Q2- 

, H _ r[l-(Q + x) 2 ] 1 

Q ~ 1 H 2 Q 2 +1, 


Wm {P + r [( Q + x) 2 - l] + x 2 - 1 } u 
{P + r l(Q + x ) 2 — 1] + x 2 — 1} u + Q 2 
Q 2 ( 1 — 2 r) — 4Qrx — 2r (x 2 — l) 

3 {P + r [(Q + x) 2 — 1}+ x 2 — l}u + 3 Q 2 ' 
2H _ 2q- 1 _ r (2 - 2(Q + x) 2 ) 1 

_ 3^2 3 3Q 2 + 3' 


(5.11a) 

(5.11b) 

(5.11c) 

(5.lid) 
(5.lie) 


Since equation (5.7e) is decoupled from the rest, we will study the stability of the reduced 
dynamical system (5.7a)-(5.7d). The scenario of mimetic F(R) gravity with arbitrary F(R) 
forms, i.e. the system of equation (5.7a)-(5.7d), admits eighteen classes of critical points (nine 
corresponding to expanding universe and nine corresponding to contracting one), where each 
class contains as many critical points as the roots of the equation M(r) = 0, with the 
exception of the curves P| which exist for the special value w m = 1 , and Pg for which 
r = —2 and M(—2) is not necessarily zero. These are presented in Table 5 along with 
their existence and stability conditions. The details of the analysis and the calculation 
of the various eigenvalues of the 4x4 perturbation matrix are presented in Appendix C. 
Additionally, for each class of critical points, using (5.11a)-(5.11e) we can calculate the values 
of various observables, such as the density parameters, the deceleration parameter and the 
total equation-of-state parameter, and we present them in Table 6 . 

The above results hold for arbitrary F(R) forms. Hence, given a specific F(R) ansatz, 
one first calculates its corresponding M(r) using (5.1) and (5.9), then he finds r* by solving 
M(r = r*) = 0, and finally one just substitutes r* in Tables 5 and 6 . 


6 Physical Implications 

In the previous sections we performed a detailed dynamical analysis for the scenario of 
mimetic F(R) gravity for exponential and power-law ansatzes, and moreover we presented 
the method for the general analysis for arbitrary F(R) forms. In this section we discuss 
on the physical features of the stable solutions, that is of the solutions that can attract the 
universe at late times, independently of the initial conditions. 

6.1 Mimetic F(R) gravity with exponential form 

The scenario of mimetic F(R) gravity with the exponential F(R) form of (3.1), exhibits 
two saddle critical points and one nonhyperbolic one, namely E]. In the latter case the 
present linear analysis is not adequate to determine its stability, and thus one needs to 
apply the center manifold method [101]. However, we mention that all the above points 
exist also in usual F(R) gravity [102], and this is explained since the extra parameter of 
mimetic gravity, namely C^, in this case is zero. Therefore, we deduce that mimetic F(R) 
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Name 

P 

Q 

r 

X 

Existence 

Stability 

pt 

0 

0 

r* 

€ 

always 

unstable (stable) for M' (r*) > 0,r* < — 1 or 

M' (r*) > 0,r* > 

pi 

0 

2e 

r* 

— € 

always 

unstable (stable) for 
—1 < w m < §, M' (r*) < 0,r* < —| or 
—1 < w m < | ,M' (r*) < 0,r* > — 1 

Pi 

9 

3 

r* 

3 

always 

stable (unstable) for 0 < w m < 1, M' (r*) <0, - l<r*<— | 

saddle otherwise 

Pi 

Pac 

2 r*£ 

V(r*-2)r*+3 

r* 

3(r*+l)« 

V(r*-2)r*+S 

’■*<-! 

numerical determination 

Pi 

0 

<?5 

r* 

>*+i(r'+2) c 

y/ r* [r* (9r*+19)+13]+4 

-1 <r*< 

or r* < —| 

unstable (stable) for 
-1 < r* < -i, M' (r*) > 0 
saddle otherwise 

Pi 

0 

2e 

3(tUm — 1) 

r* 

(ttwm-i)e 

w m<l 

stable (unstable) for 

—1 < w m <0,-1 <r* < — |(iu m + 1 ),M' (r*) < 0. 
saddle otherwise 

Pi 

0 

2 r* 

n 

r* 

3(r*+i)(«, m +i) 

n 

r* — — 2 , w m — — 1 or 

r, ^ , 5 1 - , —8(r’) 2 — 13r“—3 

-2 < r < —4, -1 < w m < A „-— 

_ 4, - m _ 6(r*) 2 +9r*+3 

or — | < r* < 0, - 1 < w m < 5 (—4 r* — 3) 

stable (unstable) for M'(r*) > 0, -1.64 < r* < —1.328, -1 < w m < 
or M'(r*) > 0, -1.328 < r* < -1, -1 < w m < 0 
saddle otherwise 

Pi 

0 

e 

r c8 

0 

1» m = i 

saddle 

PS 

0 

2 t 

-2 

0 

always 

stable (unstable) for w m > — 1 , M(— 2) > 0. 


Table 5. The real critical points and curves of critical points of the system (5.7a)-(5.7d) 
of mimetic F(R) gravity, for arbitrary F(R) asantzes. We use the notation e = ±1, where 
e = +1 corresponds to expanding universe and e = —1 to contracting one, with the stabil¬ 
ity conditions outside the parentheses corresponding to e = +1 while those inside the parenthe¬ 
ses correspond to e = —1. The symbol r* denotes the roots of the equation M(r) = 0, i.e. 
r* = M~ 1 { 0). 


Furthermore. 


we define P± c = — 2 ^ r , , 

(r*—2)r*+3 ’ 


Q 5 = 


•y/r* + l(2r* + l)e 


- v /r*[r*(9r* + 19)+13]+4’ 


fi = 


l) 2 W, 


™ “I" 6 [r* ( r * + 4) + 2] w m + (r*) 2 — 2r* + 3 and w m = ^*2 ) 3+2°y ? ) 3+ a^.. T '+ 9 j 27 “ 


4V2 

3 


48 (r ^+i36^r^ +uuHr 1 . Additionally, r c s is the parameter of curve P|. 


gravity with exponential ansatz, presents the same asymptotic behavior with standard F(R) 
gravity, and thus it does not lead to novel asymptotically late-time features. Additionally, 
note that apart from the finite critical points of Table 1, there could be stable points at 
“infinity”, which requires to apply the Poincare central projection method [103]. However, 
since this investigation lies beyond the scope of the present work, and moreover since these 
points exist also in usual F(R) gravity and thus are not new, we do not analyze them in 
more details. Finally, note that the two saddle points £2 and £3, which correspond to dark- 
energy-dominated (Ha = 1 ), accelerating (q = — 1 ) solutions, where dark energy behaves as 
cosmological constant (wde = -l), and hence they are de Sitter solutions, can be very good 
candidates for describing the inflationary phase of the cosmic evolution. 

In order to present the above behavior more transparently, we numerically evolve the 
autonomous system at hand, and in Figs. 1-3 we depict the phase-space behavior. In this 
example, the critical point £1 is the stable late-time state of the universe. 

6.2 Mimetic F(R) gravity with power-law form 

The scenario of mimetic F(R) gravity with the power-law F(R) form of (4.1), focusing in 
the more physically interesting case of expanding universe, exhibits two stable critical points, 
namely Tg~, T{g, as well as a stable curve of critical points, namely T ( 7 . We mention here that 
there are four more critical points that might be stable in a small region of the parameter space 
(namely T+, T+, T+ and T 15 ), however their exact behavior requires numerical examination. 


- 17 - 





































Name 

^m/ fn 

q 

Wtot 

n 

NWD 

NWD 

NWD 

Pi 

0 

1 

l 

3 

pe 

r 3 

0 

1 

1 

3 

pe 

-*4 

0 

— — -1 

2 r* 

- -1 

pe 

-*5 

0 

3 I Y 

2r*+l r*+ 1 1 

2 2-i 

2r*+l 3(r*+l) 1 

pe 

M3 

2 Swjyi 

1 

i 

3 

pe 

r* [r* (6ui m +8)+9v)m + 13J+3(iD m + l) 

3(ui m +l) i 

r*+Wm. + 1 

2 (r*) 2 

2r* 1 

r » 

pe 

r 8 

1 

1 

1 

3 

pe 

M? 

0 

-1 

-1 


Table 6. The real critical points and curves of critical points of the system (5.7a)-(5.7d) of mimetic 
F(R) gravity for arbitrary F(R) asantzes, and the corresponding values of the rescaled matter density 
parameter of the deceleration parameter q , and of the total equation-of-state parameter Wtot, 

calculated through (5.11a)-(5.11e). We use the notation e = ±1, where e = +1 corresponds to 
expanding universe and e = — 1 to contracting one. The symbol r* denotes the roots of the equation 
M(r) = 0, i.e. r* = M -1 (0). Furthermore, NWD stands for “Not well-defined”. 



Figure 1. The phase space of the system (3.10) of mimetic F(R) gravity with the exponential form 
(3.1), for the choice w m = 0. The point Si attracts an open set of orbits. The Figure shows the 
existence of closed orbits too. The behavior is qualitatively the same for different choices of w m . 


Point Tg is stable for n > 2, however the corresponding Q m it not well defined. Point 
T [g is a stable physical critical point, and thus it can be the late-time state of the universe. 
It corresponds to a dark-energy-dominated universe, which however is non-accelerating and 
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Figure 2. Projection of the orbits of Fig. 1 on the x-Qa plane. 



Q 

Figure 3. Projection of the orbits of Fig. 1 on the Q-fl\ plane. 


the dark energy behaves as radiation, which are not favored by observations. This point exist 
also in standard F(R) gravity [16, 104, 105] as it corresponds to C$ = 0. 

The critical points of the curve T ^ correspond to dark-energy-dominated, accelerating 
solutions, where dark energy behaves as cosmological constant wde = — 1 , and hence they 
are de Sitter solutions. They exist also in standard F(R) gravity [16, 104, 105] as they 
correspond to = 0 . 

Finally, in order to study the stability of the points T-q, Ty 2 , T-jtj and T 15 that require 
numerical investigation, we numerically evolve the autonomous system for various parameter 
choices and we depict the resulting phase-space behavior. In Fig. 4 we can see that point 
Tjt, is an attractor. Notice also the presence of heteroclinic orbits connecting the contracting 
de Sitter solution Ty 2 with the expanding one T 22 , he. corresponding to bouncing orbits 
[32, 106-109]. In Fig. 5 we observe that point Ty S is stable and thus it can attract the 
universe at late times, with the presence of bouncing solutions also visible. Lastly, in Fig. 6 
we show the stable behavior of point T[t. 

As we observe, we do find many critical points, some of which exist also in the case of 
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Figure 4. Projection of the phase space of the system (4.10) of mimetic F(R) gravity with the power- 
law F{R) form of (4-1), on the invariant set r = —n, for the choice w m = 0,n = 2. Point T+ 2 is 
the late-time attractor for the universe. Notice also the presence of heteroclinic orbits connecting the 
contracting de Sitter solution Tf 2 with the expanding one T+ 2 , i.e. corresponding to bouncing orbits. 



Figure 5. Projection of the phase space of the system (4.10) of mimetic F(R) gravity with the power- 
law F(R) form of (4-1), on the invariant set r = —n, for the choice w m = 0,n = 1.2. Point T [g is 
the late-time attractor for the universe. Notice also the presence of heteroclinic orbits connecting the 
contracting de Sitter solution with the expanding one T-j1j, i.e. corresponding to bouncing orbits. 


usual F(R) gravity, and some of which are novel and characterized by a C^-value different 
from zero. However, concerning the stable critical points, i.e. the points that can attract the 
universe at late times, we observe that they all have C$ = 0, that is they exist in usual F(R ) 
gravity too (the points that have C^ ^ 0, namely T^~, T 6 + and T 7 , are always not stable). 
This implies that, although the new features of mimetic F(R) gravity can affect the universe 
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Figure 6. Projection of the phase space of the system (4.10) of mimetic F(R) gravity with the power- 
law F(R) form of (f-1), on the invariant set r = 0, for the choice w m = 0, n = 2. Point is the 
attractor of an open set of orbits. 


evolution at early and intermediate times, that is affect the specific universe evolution, at late 
times they will not have any effect, and the universe will result at states that coincide with 
those of usual F(R). Correspondingly, the involved observables in these late-time solutions, 
do not depend on C$ either. Thus, although mimetic F(R) gravity could drive inflation in 
a different way than usual F(R) gravity, concerning the dark-energy era it cannot lead to a 
different behavior. From the dynamical system point of view this is expected, since the new 
term behaves as ~ 1/a 3 , which is known to usually lead to saddle behavior [78, 79]. Hence, 
although this term can affect the phase-space evolution, it cannot affect the stable late-time 
attractors. 

6.3 Mimetic F(R) gravity with arbitrary F(R) form 

The scenario of mimetic F(R) gravity with arbitrary F(R) forms admits eighteen classes 
of critical points (nine corresponding to expanding universe and nine corresponding to con¬ 
tracting one), where almost each class contains as many critical points as the roots of the 
equation M(r ) = 0. Amongst them, and focusing on expanding solutions, , P ( f, Pj~ 
and Pq can be stable, and thus they can attract the universe at late times. Pff and Pff 
correspond to non-accelerated universe, and thus they are not favored by observations. Pg 
is the de Sitter solution, corresponding to dark-energy dominated, accelerating universe, 
where the dark energy behaves as a cosmological constant. Additionally, point P ./~ is the 
most interesting solution, since it corresponds to dark-energy dominated, accelerating uni¬ 
verse, with dark-energy equation-of-state parameter different than — 1 , which additionally 
can have 0 < Q. m < 1 and thus it can alleviate the coincidence problem since dark energy 
and dark matter density parameters are of the same order. Finally, concerning the curve 
of points Pg that exist for w m = 1/3, physically corresponding to radiation, the fact that 
they are saddle and completely dominated by radiation energy density, may correspond to 
the radiation-dominated phase in which the universe transiently goes through its evolution, 
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before departing towards the subsequent phases. 


7 Conclusions 

In the present work we investigated the cosmological behavior of mimetic F(R) gravity. 
This scenario is the F(R ) extension of usual mimetic gravity classes, which are based on 
re-parametrizations of the metric using new, but not propagating, degrees of freedom, that 
can lead to a wider family of solutions. Indeed, in the cosmological equations one obtains a 
novel term of the form C^/a 3 , and when the new parameter goes to zero he re-obtains the 
solutions of usual F(R) gravity. In order to bypass the complexity of the involved equations 
we performed a detailed dynamical analysis, for the cases of exponential and power-law 
F(R) ansatzes, and we provided the tools to perform the analysis in the general case of 
arbitrary F(R) forms. Hence, we first extracted the critical points of the system, and then, 
for each of these solutions, we calculated various observables, such as the dark-energy and 
matter density parameters, the dark-energy and total equation-of-state parameter, and the 
deceleration parameter. 

In our analysis we found many critical points, some of which exist also in the case of 
usual F(R) gravity, and some of which are novel and characterized by a C^-value different 
from zero. However, concerning the stable critical points, i.e. the points that can attract the 
universe at late times, interestingly enough we found that they all have C$ = 0, that is they 
exist in usual F{R) gravity too. This implies that, although the new features of mimetic F{R) 
gravity can affect the universe evolution at early and intermediate times, that is affect the 
specific universe evolution, at late times they will not have any effect, and the universe will 
result at states that coincide with those of usual F(R) gravity. Correspondingly, the involved 
observables in these late-time solutions do not depend on C<j, either. Thus, although mimetic 
F(R) gravity could drive inflation in a different way than usual F(R) gravity, concerning the 
dark-energy era it cannot lead to a different behavior. From the dynamical system point of 
view this was expected, since the new term behaves as ~ 1 /a 3 , which is known to usually lead 
to saddle behavior [78, 79]. Hence, although this term can affect the phase-space evolution, 
it cannot affect the stable late-tinre attractors. 

However, we should mention that the dynamical analysis provides information for the 
background behavior only. Hence, although mimetic F(R) gravity at late times leads to 
background solutions that exist in usual F(R) gravity too, the behavior of the perturbations 
is expected to be different, since the new term contributes to the perturbations even if it does 
not contribute to the background level. Thus, it would be both necessary and interesting 
to study the effect of mimetic F(R) gravity on perturbation-related observables, such as the 
growth-index. Since this investigation lies beyond the scope of the present work, it is left for 
a future project. 

Acknowledgments 

The authors would like to thank Anupam Mazumdar, Shinichi Nojiri, Sergei D. Odintsov 
and Alex Vikman for useful comments. GL was supported by COMISION NACIONAL 
DE CIENCIAS Y TECNOLOGIA through Proyecto FONDECYT DE POSTDOCTORADO 
2014 grant 3140244 and by DI-PUCV grant 123. 730/2013. Thanks are due to all the 
members of Grupo inter-universitario de Astroffsica, Gravitacion y Cosmologia, for their 
support in a warm working environment. ENS wishes to thank Maternite Port Royal in 


22 


Paris, for the hospitality during the initial phases of this project, during the birth of his 
daughter. 

A Stability of the critical points of mimetic F(R) gravity with exponential 
form 

The scenario of mimetic F(R) gravity with the exponential form (3.1), i.e. the system (3.10), 
admits three isolated physical critical points which are presented in Table 1. In this Appendix 
we calculate the eigenvalues of the perturbation 5x5 perturbation matrix for each critical 
point. For Si the associated eigenvalues are {z\/2, — iy/2, 0,0}. Hence, it is nonhyperbolic 
with two imaginary eigenvalues, and therefore one needs to apply the center manifold analysis 
[101] , however such a study lies beyond the scope of the present work and thus we resorted to 
numerical examination (see Figures 1, 2 and 3). For £ 2,3 the eigenvalues must be obtained 
numerically, but at least one of them, with value 1.99778, is always positive. Thus, these 
two points cannot be attractors, and indeed numerical examination shows that these two de 
Sitter solutions are saddle points. 

B Stability of the critical points of mimetic F(R) gravity with power-law 
form 

The scenario of mimetic F(R) gravity with the power-law form (4.1), that is the system 
(4.10), admits 14 x 2 +1 = 29 isolated physical critical points (14 corresponding to expanding 
universe and their 14 counterparts that correspond to contracting universe, plus one more 
point without its “symmetric” counterpart) and three curves of critical points, which are 
presented in Table 3 along with their existence conditions. In this Appendix we calculate 
the eigenvalues of the 5x5 perturbation matrix for each critical point and curve of critical 
points. We use the notation e = ±1. 

For the critical points Tf the associated eigenvalues are {2e, 2e, 2e, 2e, 2e}. Thus, for 
e = +1 it is unstable, while for e = — 1 it is stable. 

For the critical points T| the eigenvalues read j ; 2e, 2e, 2e, — 2e| , and thus they 

are saddle points. 

For T| the eigenvalues are {lOe, 10c, 4e, — 2e, (4 — 6u; m )e}, and thus they are saddle 
points. 

For T| the eigenvalues write as |2e,4e, lOe, ( ' 8 "~ 1 1 °^ , (4 — 6 w m ) ej. Thus, for e = +1 

(respectively e = —1) it is a unstable (respectively stable) for w m < | and n < 1 or n > |, 
otherwise it is a saddle point. 

For T| the eigenvalues read {2e, 2e, — |e, |e, — 2w m e}, and therefore they are saddle 
points. 

For T| the eigenvalues are j —|e, — |e, , — 2w m e, 2e| , and thus they are saddle 

points. 

For Tj the eigenvalues are extracted as 

J 6(n — l)e 6ne Aie 

yjn(n + 2) + 3’ yjn{n + 2) + 3 ’ 2(n — 1) [ n{n + 2) + 3] 4 ’ 

_AA_ X^e _ 1 

2(n - 1) [n(n + 2) + 3] 4 ’ 2(n - 1) [n(n + 2) + 3] 4 J ’ 
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where A* are the three roots of the polynomial 


P( A) = A 3 y/n(n + 2) + 3 + 6 A 2 (n — 1) [n(n + 2 ) + 3 ] 4 ( 2 nw m + 1) 

— 8 A (n — 1) [n(n + 2) + 3 ] 15 ^ 2 [n 2 (32n — 9 w m — 76) + n(9w m + 51) — 9] 
—96 w m n(n — l) 2 (4n — 3)[n(8n — 13) + 3] [n(n + 2) + 3 ] 11 . 


Thus, in the general case the signs of the eigenvalues cannot be determined analytically and 
one needs to examine them numerically. For instance, for n = 2 and w m = 0 the eigenvalues 

becomes (\/lT + \/45l) , and in this case the points are sad¬ 

dle. Furthermore, since X 'f exists for n € [0.75,1.35] it follows that at least one eigenvalue 
for X 'j is positive, thus, it cannot be stable. 

For X| the eigenvalues write as 

2(n - 1 )(n — 2)e e {n [-2n(3w m + 4) + 9 w m + 13] - 3(w m + 1)} 
n* ’ n* 

[2n(7 - An) - 5]e 2(n - 2)ne [(13 - 8 n)n - 3]e 




n* 


n* 


where n* = \Jn — 1 \Jn [n(9n — 19) + 1 ]) — 4. Hence, restricting ourselves to the physical 
case — 1 < w m < 1 , we deduce that for e = +1 (respectively e = — 1 ) the point is stable 
(respectively unstable) for w m > —l and n > 2 , otherwise it is saddle. 

For Tg the eigenvalues are found to be 

' 2(3w m - l)e 2(w m + l)e 2(w m + l)e 2(3w m - 2)e 2w m e 


3{w m - 1) w m - 1 

and thus they are saddle points. 

For Tf 0 the eigenvalues are 

2(3 w m — l)e 2e(4n — 3w m — 3) 


Wr, 


- 1 


3{w m - 1) 


W r , 


- 1 


2(3 Wm 2j)d 2{Wrn 4“ 1 )^ 2jWuiC 


3(w m — 1) ’ 3(n — l)(w m — 1) ’ 3{w m — 1) 


Wm. - 1 


W m . ~ 1 


Therefore, for e = +1 (respectively e = —1) they are stable (respectively unstable) for 
w m < — 1 and \ (3w m + 3) < n < 1, otherwise they are saddle points. 

For Tf, the eigenvalues write as 

J 6 (n — l)(w m + l)e 6n(w m + l)e e {V n — l \/^2 + 3n[(2n — 3)w m — 1 ] + 3 w m + 3} 

[ \/A7 ’ \/Ai ’ 2(n— 1)\/A7 

e {y/n - 1\/Ai + 3n[(2n — 3)w m - 1] + 3w m + 3} 6 nw m e ) 


2 (n- l)v/ST 


\/AT 


where Ai = n 2 + 9(n — 1 ) 2 Wm + 6 [(n — A)n + 2 ]w m + 2n + 3 and A 2 = An 3 (3w m + 8) 2 — 
4n 2 [3u; m (18u; m + 55) +152] +3n{w m + l)(87u; m +139) — 81(tc m +1) 2 . Thus, Xq (respectively 
XQ) is stable (respectively unstable) for n = 2,w m = -1 or | < n < 2,-1 < w m < 

— 8n 2 +13n—3 

Qn' 1 — 9n+3 


For Tf 2 the eigenvalues read as 

3 e (3 y/n + y/25n^32) e [y/ny/25n^32 — 3 n) e 


0 ,- 


y/2’ 


2y[2 y/n 


2y/2, 


n 


3(w m + l)e 
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thus T+ (respectively Tu) has a 4D stable (respectively unstable) manifold for 0 < n < 
2 ,w m > — 1. However, since there exist a zero eigenvalue the points are nonhyperbolic, thus 
in order to determine their stability we need to resort to numerical examination (see Figure 

4). 


For Tf 3 the eigenvalues are extracted as 

Are 


0,-3 


n 


A26 


A36 


n + 2 ’ 2(n - l)(n + 2) 2 ’ 2(n - l)(n + 2) 2 ’ 2(n - l)(n + 2) 2 
where A* are the three roots of the polynomial 

P( A) = A 3 + 6A 2 {n — 1 )\/n(n + 2) 3//2 (ic m + 2) + 4A(n — l)n(n + 2) 3 [9(n — l)w m + 5n — 1] 
—96(n — 2)(n — l) 2 n 3 / 2 (n + 2 ) 9 ^ 2 (w m + 1). 


Hence, in the general case the signs of the eigenvalues cannot be determined analytically and 
one needs to perform a numerical investigation. For instance, for n = 2 and w m = 0 the 
eigenvalues become ^ ,0,0^*, and thus the points exhibit a 3D stable manifold. 

A complete stability analysis requires to use the center manifold theorem [101], however since 
this lies beyond the scope of the present work, we resort instead to numerical elaboration 
(see Figure 5). 

For T14 the eigenvalues write as {0, —1,1,4, 4}, and thus it is a saddle point. 

For T15 the eigenvalues write as {0, 0,0,0,0}, and thus it is non-hyperbolic. In order to 
examine its stability one needs to apply the center manifold analysis [101], however such a 
study lies beyond the scope of the present work. However, numerical elaboration allows to 
conclude that it is a local attractor (see Figure 6). 

For the curve of critical points T\§ the eigenvalues are 
jo, — 5( ^ 2 c16 ,2 Q c i 6 , — \Qd§{wm + 1) 5 — 3( ^ 2 c16 j ; where Q c \q is the parameter of the curve. Hence, 
all the points of this curve are saddle points. 

For the curve of critical points T17 the eigenvalues are 


0 , ~Qc17 


yj (75n - 32)Q 2 17 - 25n 

\A^\/3Qc17 ~~ 1 


- 3 


:Qci7 


\/ (75n - 32)Q 2 17 - 25n 

V^V^Qc17 ~ 1 


- 3 


1 3Qcn{w>m 4“ l)j 3Q c 17 


(B.l) 


where Qd7 is the parameter of the curve. Therefore, the corresponding points are stable for 
the combinations: 


n < 0, —1 < w m <1,5 


n 1 

75n - 32 < ^ c17 < ^3 


2 ^ 1 
or0 <n<-,— 1 < w m < 1 , Qci 7 > ^ 

2 1 ^ I n 

orn>-,-l<u; m <l,- 7 =<Q cl 7<^^-^. 


(B.2) 


For Tl 8 the eigenvalues read as | — |e, — |e, — —{w m + l)e, —e j. Thus, for e = — 1 

(respectively e = +1) it is unstable (respectively stable) for w m > —1 and 0 < n < 1, 
otherwise it is a saddle point. 
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C Stability of the critical points of mimetic F(R) gravity with arbitrary 
F(R) forms 

The scenario of mimetic P(P) gravity with arbitrary potentials, i.e. the system of equa¬ 
tion (5.7a)-(5.7d) admits eighteen classes of critical points (nine corresponding to expanding 
universe (e = +1) and nine corresponding to contracting universe (e = —1)), where each 
class contains as may critical points as the roots of the equation M(r = r*) = 0, with the 
exception of the curves P| which exist for the special value w m = and Pg for which r = —2 
and M{— 2) is not necessarily zero. These are presented in Table 5 along with their existence 
conditions. In this Appendix we calculate the eigenvalues of the 4x4 perturbation matrix 
for each critical point and curve of critical points. 

For the critical points P( the associated eigenvalues are |2e, 2e, 2 ^. r «+i^ 6 ', 2 eM' (r*)j. 

Thus, for e = +1 (respectively e = — 1) they are unstable (respectively stable) for M' (r*) > 
0, r* < — 1 or M' (r*) > 0,r*> — A, otherwise they are saddle points. 

For the critical points P| the eigenvalues read |4e, 2 ^»_ t f 1 >>e , —2(3 w m — 2)e, —2 eM' (r*) j. 
Therefore, for e = +1 (respectively e = — 1) they are unstable (respectively stable) for 
— 1 < w m < |< 0, r* < —| or —1 < w m < | ,M'(r*) < 0, r* > —1, otherwise they 
are saddle points. 

For the critical points P| the eigenvalues write as j — 2 3 (v*+ij 6 ’ —2 w m e, | eM' (r*)j. 

Therefore, for e = +1 (respectively e = — 1) they are stable (respectively unstable) for 
0 < w m < 1, M' (r*) <0, — 1 < r* < — |, otherwise they are saddle points. 

For the critical points the eigenvalues are the roots of the polynomial 


P(A) = A 3 A a (r*) 3 - (r*)' z + r* + 3 




+ 3A 2 [(r *) 3 - (P ) 2 + r* + 3j ( 2 r*w m - 1) 

—2A 3 A {r* [r* (32 r* + 9 w m + 76) + 9 w m + 51] + 9} 

-384 (r *) 4 w m - 912 (r *) 3 w m - 612 ( r*) 2 w m - 108 r*w m (w m + 1), 

where A 3 = \/(r *) 2 — 2 r* + 3, and the fourth eigenvalue is _ 6 h +i)M (r ) ^ jj ence) j n 
V y (r*—2)r*+3 

general case the signs of the eigenvalues cannot be determined analytically and one needs to 

perform a numerical investigation. 

For Pp (r*) the eigenvalues are the roots of the polynomial 


P(A) = A 3 A a (r*) 3 - (r*) z + r* + 3 




— 3A 2 (r*) 3 - (r *) 2 + r* + 3 (2r*w m - 1) 

—2A 3 A {r* [r* (32 r* + 9 w m + 76) + 9 w m + 51] + 9} 

+384 (r *) 4 w m + 912 (r *) 3 w m + 612 ( r *) 2 w m + 108 r*w m (w m + 1), 

and the fourth eigenvalue is 6 h +^) M ( r ) _ xi 1USj ; n the general case the signs of the eigenvalues 

y (r*-2)r*+3 

cannot be determined analytically and one needs to perform a numerical investigation. 

For P| (r*) the eigenvalues are 


(2r* + 1) (4 r* + 5) e 


( 2 r* + 1 ) (4r* + 5)e 


\Jr* + 1 \Jr* (r* (9 r* + 19) + 13) + 4 y/r* + 1 y/r* (r* (9r* + 19) + 13) + 4 
(r* (8r* + 13) + 3) e 2 y/r* + 1 (r* + 2) eM' (r*) 


y/r* + 1 y/r* (r* (9 r* + 19) + 13) + 4 y/r* (r* (9 r* + 19) + 13) + 4 


}• 
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Thus, P 5 + (respectively P 5 _ ) is unstable (respectively stable) for — 1 < r* < — ^ , M 1 (r*) > 0 
or r* < — |, otherwise they are saddle points. 

For the critical points P| the eigenvalues write as 

f 2w m e 2e (4r* + 3 w m + 3) 2(3w m — 2)e 2(3w m - l)eM' (r*) 1 

I w rn, - 1 ' 3 (r* + 1) (ic m - 1) ’ 3{w m - 1) ’ 3(w m - 1) J ' 

Therefore, for e = +1 (respectively e = — 1) they are stable (respectively unstable) for 
— 1 < w m <0,-1 < r* < — f (w m + 1), M' (r*) < 0, otherwise they are saddle points. 

For the critical points P| the eigenvalues write as 

f 6 r*w m e -y/A[} + 6 r*w m e + 3(w m + l)e + 6 r*w m e + 3{w m + l)e 

I VSI ’ 2 VAI ’ 2y/Kl 

6 (r* + 1) (w m + l)eM' (r*)\ 

VAl J ’ 

where 


A 4 = 9 (r* + if w 2 m + 6 [r* (r* + 4) + 2] w m + ( r*f - 2r* + 3 
A 5 = (r* + I)' 1 {r* {4r* [r*(3w m + 8) 2 + 3w m (18w m + 55) + 152] 

+3 (w m + l)(87w m + 139)} + 81(«7 m + l) 2 } . (C.l) 


Thus, Pf (respectively P 7 ) is stable (respectively unstable) for M'(r*) > 0, —1.64 < r* ^ 

where W~ = -32(r*) 3 -110(r«) 2 -113r«-27 4^2 H 

J, wnere tt> m 3[4(r*) 3 +24(r*) 2 +29r*+9] 3 ^ 


— 1.328, —1 < u; m < w r 


48(r*) 5 +136(r*) 4 +115(r*) 3 +25(r* ) 2 


[4(r*) 3 +24(r*) 2 +29r*+9] 


or M'(r*) > 0, —1.328 ;$ r* < —1, —1 < w m < 0. It is a saddle otherwise. 

For the critical points P| the eigenvalues write as {4e, — e, e, 0}. Thus, they are saddle. 
For the critical points P| the eigenvalues are 


3(w m + l)e 

V2 


3e 

V2’ 


(y9-8M(—2) + 3J 

le 1 

(y9 - 8 M(—2) - 3J 

"1 

2 y /2 

5 

2V2 

I 


Thus, Pf (respectively P 9 ) is stable (respectively unstable) for w m > —1,M(—2) > 0. 
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